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Abstract 

We analyze the Ising model on a random surface with a boundary magnetic field 
using matrix model techniques. We are able to exactly calculate the disk amplitude, 
boundary magnetization and bulk magnetization in the presence of a boundary field. 
The results of these calculations can be interpreted in terms of renormalization group 
flow induced by the boundary operator. In the continuum limit this RG flow corre- 
sponds to the flow from non-conformal to conformal boundary conditions which has 
recently been studied in flat space theories. 
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1 Introduction 



Simple statistical mechanical lattice models, such as the Ising model, have been used for 
many years to gain insight into the behavior of a wide range of physical systems. The great 
utility of such simple models arises from the fact that in many cases they are exactly solvable. 
It is a remarkable fact that some questions which seem analytically intractable for models 
on a regular lattice yield exact solutions when the underlying lattice is itself taken to be a 
random element of a larger statistical ensemble. An example of such a situation is the Ising 
model with a bulk magnetic field. Although this model has not been analytically solved on 
a fixed lattice, it is possible to exactly compute the partition function and magnetization of 
the model on a random lattice 

In this paper, we consider the Ising model on a random lattice in the presence of a 
magnetic field on the boundary, rather than in the bulk. Again, this corresponds to a 
problem which does not seem to be analytically solvable on a fixed lattice. Summing over 
random lattices, however, we find that the partition function and magnetizations can be 
calculated exactly. 

The Ising model on a random two-dimensional lattice can be described in terms of a 
matrix model. A great deal of technology has been developed to deal with matrix models, 
primarily as a tool for studying string theory. The techniques we use here were derived in 
an earlier pair of papers [|, ^, and are related to methods described in ||. Some of the 
results which we describe here appeared in a previous letter 0. 

The Ising model on a regular lattice with a boundary magnetic field was studied many 
years ago 0. This model has been of renewed interest recently because in the continuum 
limit, where the lattice spacing is taken to zero, it gives a simple example of a two-dimensional 
field theory which is conformal in the bulk but which has boundary conditions breaking 
conformal invariance P, ^, [10[- The only boundary conditions for the Ising model which 
preserve conformal invariance are free boundary conditions (where the boundary field h 
vanishes), and fixed spin boundary conditions (where h = ±oo). Putting an arbitrary field 
h on the boundary generates a renormalization group (RG) flow from the free boundary 



condition to the fixed boundary condition |T^. 

For a matrix model corresponding to fields on a random surface, a continuum limit can 
also be taken. In this limit, the theory describes conformal matter fields coupled to 2D 
quantum gravity. In this paper we take the continuum limit of the disk partition function 
and magnetizations and consider the implications of our results for the resulting theory of 
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c = 1/2 matter coupled to gravity. In particular, we find that the results are in accord with 
the hypothesis that the RG flow which has been understood in flat space is present in an 
appropriate form in the theory with gravity. 

One provocative feature of our results is that as the boundary magnetic field is increased, 
except for a jump discontinuity when the field becomes nonzero, the expectation value of 
the magnetization of a randomly chosen spin in the bulk decreases. This counterintuitive 
result may be explained in terms of the effects of the matter fields on the geometry — 
roughly speaking, the increase in magnetic field produces a long "throat" which separates 
the boundary from the bulk and which increases the average distance of a bulk spin from the 
boundary, effectively decreasing the bulk magnetization. However, this result is also found 
to be a finite volume effect which may depend upon the precise choice of how the random 
lattice ensemble is defined. 

Another context in which this work may be relevant is the current discussion of D-branes 
in string theory (see for example ||13|). Just as there are two conformally invariant boundary 
conditions for the Ising theory, a conformal field theory of a single bosonic field can have 
two conformally invariant boundary conditions: Neumann and Dirichlet. Considering the 
continuum limit of the Ising model as a single free fermion, free and fixed Ising boundary con- 
ditions are related through supersymmetry to Neumann and Dirichlet boundary conditions 
on a bosonic field. The fact that RG flow behaves similarly in flat space and in the presence 
of a fluctuating metric suggests that perhaps Dirichlet boundary conditions in superstring 
theory naturally arise as an RG limit of a non-conformal boundary term. 

In Section 2 of this paper we describe the discrete model we will use for the Ising model on 
a random surface. In Section 3 we apply the methods of |@] to this model, explicitly deriving 
a quartic equation satisfied by the disk amplitude and locating the associated critical point. 
In Section 4 we take the continuum limit of the disk amplitude. In Section 5 we discuss the 
general formalism we will use for calculating magnetizations, and we apply this in Sections 
6 and 7 to compute the boundary and bulk magnetizations on the disk. In Section 8 we 
discuss implications of our results for renormalization group flow, duality, and the effects of 
gravity on the behavior of the matter theory. 

2 The model 

A discretized theory of c = 1/2 matter coupled to 2D quantum gravity is described by the 
Ising model on a randomly triangulated surface. At the center of each (equilateral) triangle 
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on the surface lives a single Ising spin, coupled to its nearest neighbors. This theory can be 
described by a matrix model [|I| with partition function 

Z{g, c) = JdU BV exp {-NS{U, V)) , (1) 

where the action is given by 

S = -Tr + V^) - cTr UV - ^Tr (f/^ + V^) . (2) 
2 3 

In these expressions, U and V are N x N hermitian matrices representing up and down spins 
respectively, g is a coupling constant corresponding to a Boltzmann weight for each triangle 
on the surface, and c describes the coupling between Ising spins. (This c is unrelated to the 
central charge.) The partition function can be expanded in a power series in g and and 
the coefficient of g^N"^'"^^ is then given by a summation of the Ising partition function over 
all triangulations of a genus h Riemann surface by k triangles. 

We are interested in amplitudes corresponding to various boundary conditions on spins 
living on a genus zero surface with boundary. For example, the disk amplitude for a configu- 
ration of plus and minus spins on the boundary, represented by an ordered string w{U, V) of 
the matrices U and V , is given by the large- limit of the matrix model expectation value, 

hm l(Tr^^/,r)) . (3) 

(Since we will always be interested in the large- limit, henceforth we will suppress the 
explicit in expectation values.) Again, can be expanded in a power series in g, with 
the coefficient of g'' giving the sum over all disk triangulations by k triangles with a boundary 
having a fixed length and spin configuration w{+, —). 

There are two types of conformally invariant boundary condition in the Ising theory: 
"fixed" (corresponding to all boundary spins aligned) and "free" (corresponding to an equally 
weighted sum of all possible boundary spin configurations). For each of these we can define a 
generating function which encodes the amplitudes. Fixed boundary conditions are described 
by 

oo 

^fixed(n) = ^(Trf/'^)M\ (4) 
while free boundary conditions are described by 

oo 

0free(x) = 5:(Tr(t/ + V)^)a:^ (5) 

k=0 
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The Ising model in this set of variables is symmetric under interchange of U and V (corre- 
sponding to the symmetry under spin reversal). The generating function for fixed boundary 
conditions can therefore also be described by summing over amplitudes with all V^s on the 
boundary. 

An alternative coupling of discretized c = 1/2 matter to 2D quantum gravity is obtained 
by considering the dual Ising model on a random surface [Q. The partition function of 
the dual model is defined once again as a sum over surfaces with Ising spins at the face of 
each plaquette, where now the plaquettes may be polygons with arbitrary numbers of sides, 
but the coordination number at each vertex is constrained to be equal to three. Such a 
configuration is equivalent through duality to a triangulation in the original theory, but with 
spins located at vertices rather than on faces. The dual model may also be described as a 
theory of two matrices X and Y, with action 

S = ^^^Tr X^ + ^^^Tr Y^ - ^Tr (X^ + SXF^) . (6) 

There are once again two types of conformally invariant boundary conditions, fixed and free. 
In the dual model the matrix variables X and Y do not refer to the state of individual spins, 
but rather to the relative state of two spins across an edge; X denotes an edge separating 
two equal spins, while Y denotes a boundary between two opposite spins. The generating 
function for fixed boundary conditions is therefore 

oo 

0fixed(x) = E(TrX'=)a;\ (7) 

fc=0 

while free boundary conditions are described by 

oo 

0free(w) = E(Tr(X + F)'=)M^ (8) 

Note that the dual model is symmetric under Y —Y\ it is therefore the generating function 
for free boundary conditions which has an alternative description in these variables, as a sum 
over (Tr (X — y)'^). (Amplitudes with an odd number of boundary y's vanish identically.) 

Although the original and dual formulations of the Ising model represent different cou- 
plings to 2D gravity, there is a simple transformation that relates the original action (H) to 
the dual action (H): 

X ^^iJJ^V) 
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1^2 
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As a consequence, the partition functions for the two models on surfaces with no boundaries 
are identical. This does not, however, guarantee that all correlation functions in the two 
theories agree, since the transformation has a nontrivial action on the states and operators 
of the theory. For example, the disorder operator Y in the dual theory is taken into the spin 
operator U—V in the original theory. Similarly, the role of free and fixed boundary conditions 
is interchanged; we have 

0frcc(a;) = 0fixed(V2a;) . 

The (Kramers-Wannier, or T-) duality of the model relates a specified state (fixed, free, 
or with additional operator insertions) in the original form of the model to the same state 
in the dual version. Although in the discrete version of the theory, this duality symmetry 
is explicitly violated by finite volume effects, it seems likely that the duality symmetry is 
restored for all correlation functions in the continuum limit. Evidence for duality in the 



continuum limit at genus zero was presented in |15, p[, and the issue of higher genus was 



explored in [16, 17, 18 



Our interest in this paper is in non-conformally-invariant boundary conditions, which 
may be thought of as arising from the introduction of boundary fields or couplings. We 
therefore consider the original model in a new set of matrix variables Q, R, defined by 



Q = e^U + e-^V , 
R = e^U - e-^V . 

Substituting these into the original action (||), we obtain 

S = \Tt [(cosh(2/i) - c)Q2 + (cosh(2/;,) + c)R^ - 2 smh{2h)QR] 
-^Tr [cosh(3/i)(Q3 + SQR^) _ smh{3h){3Q^R + R^)] . 



(11) 



(12) 



In this set of variables we can calculate the generating function for disk amplitudes with all 
Q's on the boundary, 

oo oo 

0(g,/z) = E(Tr {e'^U + e-'V)V = E(Tr Q')?' • (13) 

k=0 fc=0 

This generating function describes boundary conditions which interpolate between fixed and 
free. The parameter h can be thought of as a boundary magnetic field applied to otherwise 
free boundary conditions; for h = we recover (^free in the original model, while for h = ±oo 
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the boundary spins are all driven to one value and we recover ^gxed- Note that although the 
h ±00 limit appears to be singular, this is just an issue of normalization, which can be 
absorbed by a constant rescaling of q. Indeed, we shall see that although the critical value 
of q goes to zero as h —>■ ±00 with the chosen normalization, all quantities of interest (such 
as for example) are well behaved in these limits. 

3 Loop equations for correlation functions 

The process of calculating the generating function (f){q) was described and essentially carried 
out in 0, without the explicit solution being written down. Here we will review the basics 
of that procedure, and examine the solution in detail. (The discussion in was framed in 
terms of non-commuting variables; for our purposes here we may skip directly to functions 
of a single variable.) 

We begin by defining an additional set of generating functions 4>w{q,r){q)- These functions 
describe disks whose boundaries include a fixed string of matrices w{Q, R), plus any number 
of additional Q's: 

00 

'^K5,r)(g) = E(Tr wiQ,R)Q')q'' . (14) 

k=0 

For example, we have 



Jrqr 



(g) = Y.i'^T RQRQ')q' (15) 

fc=0 

Prqr ~l~ QPrqrq ~l~ •? Prqrqq ~l~ •? Prqrqqq ~l~ . . . . (-^6) 

Recall that Pw{q,r) is the amplitude for a disk with boundary specified by the string w{q, r); 
thus Pg corresponds to a single boundary edge labelled q, Pr corresponds to a single boundary 
edge labelled r, and Po = I corresponds to no boundaries. 

We will also introduce a derivative operator Dg which acts on power series in q. Its effect 
is to annihilate terms independent of q, and to remove one power of q from all other terms: 

^.^'-C- till: (1^) 

The action of on a generating function is to remove any constant term and divide the 
remainder by q. For example, 

D,0 = q-\<p-l) , 

= , (18) 

D<?0r = q~^{4>r -Pr) , 
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Q Q 



Figure 1: Decomposition of 0^- Removal of the external edge marked R on the left hand 
side leads to one of the possibilities shown on the right hand side, depending on whether 
that edge was connected to another exterior edge or an interior triangle. 

and so on. 

We can now derive a set of loop equations which relate these functions to each other, by 
considering the possible outcomes of removing a marked edge on the boundary. For example, 
let us examine the effect of removing the edge marked R from the disks represented by (priQ)- 
Since (pr represents a sum over various triangulated geometries of the disk, we can consider 
the effect of removing R from each of the terms separately. For each term, there are two basic 
possibilities: the edge might be identified with another edge elsewhere on the boundary, or 
it might belong to a triangle. In the first case, removing the marked edge and the one it 
was connected to leaves two disconnected triangulations, both with all Q's on the boundary 
and therefore representing (f){q). In the second case, removing the triangle reveals two new 
boundary edges, which may be marked with two Q's, two R's, or one Q and one R; these 
alternatives relate the initial generating function to 0^0, (prr, and Dq(j)r- This decomposition 
of (pr is shown schematically in Fig. 1. The amount which each term contributes to (pr can 
be derived from the action ([1^), as detailed in 0. 

By this procedure we can derive a closed system of eight independent equations in the 



8 



quantities 



fri frri frqri frrri frqqri frqrri frrrrj 



1 + Pq^(f)^ + aq(f)rr + aDg0 + 2d(f)r 

2 + b(f)rr + + 2eDg(f)r 



Bg(j)r = 7Pg0 + + l3q(t)Dq(j)r + (i^^^^r + O^rgrr + aDg</>r + dViq 



(19) 



D|}0rr = 7Pr0 + 70r + Pq4>4>rr + d(j)rqrr + (Kprrrr + C^D^^j-r + dDq(j)rrr 

(j)rr = Otcj) + '^q4>4'r + C^rgr + fe^rrr + &Dg0r + cDg^rr- 

0r-gr = + 70r + '^q(l)Dq4>r + e0rqgr + fe0-rgrr + &Dg0r + eDq4>rqr 

(prrr = O.Pr<P + Ot<Pr + iq<P<Prr + G<Prqrr + b(j)rrrr + &Dg0„ + cDqlprTr ■ 

(In we listed ten equations in these variables, but the last two were not linearly 
independent from the first eight.) Here we have defined the new variables 

a = Y3^[cosh(2/i) — c] 
/3 = -^[cosh{2h) + c] 



l-c 
2 



7 = T^sinh(2/i) 



a 



2(Y^[cosh/i + ccosh(3/t)] (20) 
b = — 2[i-c^) [sinh h — c sinh(3/j,)] 
d = — 2(i^e2) [sinh h + c sinh(3/j.)] 
e = 2(i-c2) [cosh /i — c cosh(3/t)] . 

The set of equations (|19]) is completely algebraic, since we can replace the derivative operators 
with algebraic functions as in (pISl). 



Note that we can cut down somewhat on the number of undetermined variables by looking 
at equations (p!9| ) order by order. These give the following relations: 

Pq = aPrr + CiPqq + 2(ipgr 
Pqq ~l~ (^Pqrr ~l~ ^Pqqq ~l~ '^dpqq^ 

Pr = bprr + bpqq + 2epqr 

[21) 

Pgr = 7 + bpqrr + bpgqq + 2epqqr 
Pqr = 7 + ^dpqrr + aPrrr H" ^Pqqr 
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Furthermore, each of these amphtudes may be expressed as a sum of amphtudes in the 
original Ising model; for example, Pq = e^pu + e~^Pv = 2(cosh/i)p„. We use these relations 
to eliminate everything but pi = pu and ps = Puuu- 

It is now a straightforward but tedious exercise to solve the system ([l9| ) by deriving a 
single quartic equation for as a function of g, c, pi and p^. The quartic is given by 

/40' + /30' + /20' + /l0 + /O = , (22) 

with 

/4 = 2g^[l - cos\i{&h)]q\ 

/s = -Ag^ cosh(3/i)g^ 

+V[pQgl^(4/^) ^ cosh(2/i) + c(3 - cosh(6/i))]g^ 

+2 c/[c2 cosh (9/i) + ccosh(7/i) - (1 + Sc^) cosh(3/i) - (1 + 5c) cosh(/i)]g^ , 

/2 = -C/V + V[cosh(/i) -2ccosh(3/i)]g=^ 

+^2[4c2 cosh(6/i) + 12ccosh(4/;,) - 2(1 - 5c) cosh(2/i) - (3 - 2?,c^)]q^ 
+2g[-2c^ cosh(7/i) - 2c(l + 2c) cosh(5/i) - (5c + 13c=^) cosh(3/i) 

+ (1 -4c- 19c2)cosh(/i)]g5 
+2[c3 cosh(8/i) + c(2c + 2c^ + g^) cosh(6/i) 

+ (c + 2c^ + 5c^ - 5-^ + 2pig^) cosh(4/;,) 

+ (2c + 4c2 + 6c3 + ^2 - 2pi^3) cosh(2/i) + c(l + 4c + 2c^ - g^)]q^ , 

/i = -2c^^ + 2c^3[5cosh(/i) + ccosh(3/i)]g 

-2c^2[2ccosh(4/i) + 4(1 + c) cosh(2/i) + (5 + 3c2)]g2 
+2c/[c2(l + c) cosh(5/i) + c(l + 6c - 5c=^ + 3g^) cosh(3/i) 

+ (6c + 2c^ - Ac^ - fi-^ + 2pig^) cosh{h)]q^ 
+2[-c2(c - c^ + 5(2) cosh{Qh) + 2c(-c + c^ - 2^(2 + p^g^) cosh(4/i) 

+ (-c - 2c^ + + 2c^ + - 2pi5(3 - ^cpxg^) cosh(2/i) 

+ (-c + c^ + - h(?g^ - 2pi5(2)]g'^ 
+2(7[c2(l — pxg^ cosh(7/i) + c(l — p\g + 2cpig) cosh(5/i) 

+3c(c2 + pi5() cosh(3/i) 

+ (-1 + 4c2 + 2^15 + 2cpig + c^pig) cos\i{h)]q^ , 
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/o = 2c^^ + 2cg^[-ccosh{3h) + (-5 + 2pig) cosh(/i)]g 

+g^[2c^{2 - Pig) cosh(4/i) + 2c(4 + 4c - Spigf - Scpigf) cosh(2/i) 

+ (10c - 6c^ + Apig - lAcpig - 3g^ - 'ip3g^)]q'^ 
+2g[c^{—l — c + Pig + cpig) cosh(5/;,) 

+c(— 1 — 6c + 5c'^ — 3pi(7 + IScpi^f + ?>g^ + 4^3(7'^) cosh(3/;.) 
+ (-6c - 2c2 + 4c2 - Apig + IScpigf + 2c^pig 
+c'pig + S^f^ + Ap^g^) cosh(/i)]g^ 
+ [2c^(c — c'^ + pi^f — ?>cpig — g^ — Psg^) cosh(6/i) 

+2c(2c — 2c^ + 2pig — Gcpig — g"^ — Pig^ — 2p^g^) cosh(4/i) 
+2(c + 2c^ — c^ — 2c^ + Pig — cpig — Qc^pig — g^ — 2cg'^ + pig^ 

+cpig^ - psg^ - 2cp^g^ - pjg'^) cosh(2/i) 
+ (2c — 2c^ + 2pig — Gcpig + 2c^pig — Gc^pig — g^ — c^g"^ 
-2pig^ - 2p3g^ - 2c^p^g^ + 2plg^)]q' . 



(23) 



The analytic solutions to such an expression are of course rather unwieldy, but fortunately 
they are also of little interest to us. Instead, we are interested in the expansion of around 
the critical point of the model, which encodes the continuum limit of the theory. (See for 
a discussion of the extraction of the continuum limit.) The critical values of the quantities 
c, g, pi and p^ are well known [|l], 



(24) 



3-^/2^ (-l + 2v^)\ (25) 



-v^), (26) 

and 
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^cP3c = Y^(-699 + 40-7i) . (27) 

We would now like to find the critical value of q for any given boundary field h. The 
critical point is defined as the radius of convergence of the power series expansion for the 
generating function, interpreted physically as the point where boundaries with an infinite 
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number of segments begin to dominate. The generating functions for fixed boundary condi- 
tions, 0fixcd(w) in the original model and 0fixed(a;) in the dual model, obey cubic equations, 
and the critical point is simply the value of x or m for which the cubic has repeated roots. 
For the quartic this story is slightly more complicated. 

Fig. 2 shows a numerically generated plot of the real part of the solutions to the quartic, 
as a function of g, for exp(/i) = 5/3 and the other parameters at their critical values. Only 
three distinct curves are visible, as two of the solutions have identical real parts. At three 
points on the graph the curves intersect, representing multiple roots; as q increases, there 
is first a triple root, then a double root, and another triple root. To decide which of these 
represents the actual critical point, we follow the behavior of the physical solution as q is 
increased from zero along the real axis, and look for a branch point (which will indicate the 
radius of convergence). From the definition of the generating function we know that the 
physical solution is the one which equals one at g = 0. As g is increased, the first triple root 
does not represent a branch point for this solution, and is therefore not the critical point. A 
similar situation holds for the double root; even though the physical solution is one of the 
double roots, one can verify numerically that circling the double root in the complex q plane 
does not take you to a distinct Riemann sheet, and hence this value is not a branch point. 
The critical value of q is therefore at the second of the two triple roots, as indicated by the 
point (c) in the figure. 

To discover an analytic expression for the critical value gc, we note that at this point the 
quartic can be factored into the form 

/4(0-0«)(0-0(3)f = 0, (28) 

where (f)^^^ is the triple root and </)*^^-* is the single root. Setting the coefficients of ( ^81) equal 
to those of ( p2D yields a set of equations from which we can eliminate (p^^^ and 0'^'^^ to obtain 



a single octic equation for g^ as a function of h. Happily, this octic factors into the product 
of two quadratics (one of which is cubed). It is most conveniently written in terms of qd Qc- 



[729 + {zi cosh 3/i + Z2 cosh h) [q^/ g^) + (-23 cosh Ah + 2z^ cosh 2h + z^ (g^/ Qc 

X [729 + {zi cosh 3/i + Z2 cosh h){qc/gc) + (-25 cosh4/i + 22:5 cosh 2h + ZQ){qc/ gc)^] , 

(29) 
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0.16 



Figure 2: The four roots of the quartic, as a function of q. The dashed hnc represents two 
roots which are complex conjugates of each other. Moving from left to right, we find the 
triple root "a", the double root "b", and the triple root "c" which is the critical point. The 
physical branch is the one that goes to 1 as g goes to zero. 

where 

zi = 54(1-2^7) 

Z2 = -1458 

z, = 6(10 + 7^7) ^^^^ 

Z4 = 2(262 -llv^) 

Z5 = 18(-16 + 5v^) 

ze = 2(784 - SSV?) . 

The critical point is the solution obtained by setting the cubed quadratic to zero; one root 
of the quadratic will be the critical point for h > 0, and the second for h < 0. For h > the 
critical value is 

_ yc(l + 2v/7)e^^ 

^ ~ 1 + (-1 + V7)e'^ + (2 + V7)e^f^ ' ^ ^ 



whereas for /i < 0, 



^^^^ " 1 + (-1 + V7)e-^^ + (2 + V7)e-''^ ' ^ ^ 



so that although qdh) is continuous at h — 0, it is non-analytic there. 
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The corresponding critical value of (j) is 

(3e2^ - 1)(1 + (-1 + V7)e2'^ + (2 + ^)e^^) 



(33) 



Note that as expected, Qc a.s h oo, but 0c is finite and non-zero for all h. 



4 Continuum limit and disk amplitude 

To extract information about the theory in the continuum limit, we analyze the behavior of 
the discrete model in the vicinity of the critical point. The coupling constant c is set to the 
critical value Cc given by (P¥); deviations from this value would correspond to introducing a 
mass for the Majorana fermion in the continuum limit of the Ising model, which we do not 
examine in this paper. We then trade the independent variables g and q for new variables t 
and z, defined by 

9 = 9ce~ 
q = qcC -- 

where e is a small parameter indicating the distance from the critical point. The expansions 
of pi and p3 in powers of e are then given by 



(34) 



3-V7 



1 - 



22 + lOv^ 
9 



^ 5V3(55 + 25v/7) ^3/3^^/3 



^ 5^/3(11 + 5x/7) ^,o/3^,/3 



216 



36 



(35) 



and 



-699 + 40-73/2 4(121-5-73/2 



9P3 



100 
9 - 51/3 



25 



+ — - — e^'-'r'" + — - — e'^'-'t"'^ + 0{e^r) . (36) 

These expansions may be substituted into the quartic (0), which may then be solved 
for as a sum of increasing powers of e. We obtain 

3 + 2(-2 + v^)e2^ + (1 + 2V7)e^'' 



+ 



5 - 27/3a(/i) 



,4/3$ + 



lOe^^ ^ 
2^/3 (1 + 4e2'^ + e 



eZ 



4h 



15a{h) (l-e^^^) 



g5/3 2 (2T + Z2) $ 



$2 - (4r)4/3 



(37) 
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where for convenience we have rescaled the variables to 



T = 5t (38) 

and 

a{n) 



(39) 



and the function $(Z, T) is given by 



4/3 , r— —\i/3 



<I>(Z, T) = i^Z + VZ^ - ATj +[Z - _ 4Tj . (49) 
For /i > 0, 

a(/i) = ^ ^ . (41) 

^ ' 1 + (-1 + v^)e2'^ + (2 + v^)e4'* ^ ' 

but this function is discontinuous at /i = 0. If we take the continuum limit in (|22|) after /i is 

set to zero (cf. the calculation in H]), we find that cj) is given by (|37D with 



The universal part of is the first non- analytic term, and appears at order e^^^. By 
virtue of the discontinuity in a{h), both the numerical coefficient of the universal term, and 
the amplitude of as a function of z, are discontinuous. 

The universal part can be converted into the asymptotic form of the disk amplitude 
0(/, a) for fixed boundary length I and disk area a. These forms of the amplitude are related 
through a Laplace transform 

3 . ^y^^^^^ {z/a{h), 5t) = JdlJda e-^^-*"<^(/, a) . (43) 

Inverting the Laplace transform, we have 

0(/, a) = ^^(a(/.)/) V^(a/5)-^/^e-^("(^»V^ = ^^^L^/^A-y^e-^'/^ , (44) 

with the rescalings 

L = a(h)l, A = - . (45) 
5 

Up to an irrelevant multiplicative constant, this is precisely the form of the disk amplitude 
when the boundary conditions are conformal |jl9|, ^ |^ (i.e., with h = ot h = ±00); 
however, the boundary length / is rescaled by the factor a{h) which depends discontinuously 
on the boundary magnetic field. Note that this amplitude includes an extra factor of / 
corresponding to a marked point on the boundary. 
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5 Expectation values 



Now that we have defined the continuum hmit of the model and computed the disk amphtude, 
we would like to calculate a number of correlation functions related to the boundary and 
bulk magnetizations in the theory. In this section we describe the formalism necessary to 
perform such calculations efficiently. 

An example of the type of calculation we need to perform is the limiting value of the 
boundary magnetization on a disk with a large number of triangles and boundary segments. 
The boundary magnetization for a spin on the boundary of a disk with k boundary edges 
and n triangles is given by 

_ (Tr {e^U - e-''V){e''U + e~^V)^~^)n 
(^)".^ - (Tr [ei^U + e-hV)% ' ^ ' 

where by ()„ we indicate a sum over triangulations restricted to geometries with n spins (the 
coefficient of in an expansion in g). The quantity (m) is defined to be the large n and k 
hmit of (H). 

More generally, we define the expectation value of an operator in a specified state ip as 

{A)= hm i^M:^, (47) 

where {4')n,k is taken to mean the sum over all triangulations, with appropriate weights, 
with n triangles and k boundary edges. {A4')n,k is the same quantity, but with the weights 
adjusted by the operator A. 

For the cases we are considering in this paper, ip is a. sum over all triangulations, with 
weights determined by the boundary magnetic field. Thus, for zero field, ip is simply a 
sum which is equally weighted for all boundary configurations (free boundary conditions), 
whereas for infinite (positive) h, the only boundary configurations with non-zero weight are 
those with all spins pointing up (fixed boundary conditions). When A represents a boundary 
spin operator, for example, then for each configuration the boundary spin is evaluated at a 
particular site, and the weight acquires a ±1 depending on whether that spin is up or down. 
When ^ is a bulk spin operator, the spin is evaluated at a site in the bulk. 

The limits in (^7]) can be understood in terms of the asymptotic behavior of {Aip)n.k and 
{4')n,k- For large n, k these functions scale asymptotically as 

Wn,kr-g;''q;'f{n,k) (48) 
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and 

{A^p)n,k^g:y-''g{n,k) . (49) 

Thus 

(^). ,„„ <p]^ (50) 

For large n and k, it is appropriate to replace the number of triangles and boundary edges 
by the area and length variables a = e^n and / = ek, so that (A) will appear as a function 
of a, I. 

The continuum limit of an operator expression of this type can be easily determined by 
taking the continuum limits of the quantities {iJj) and {Aijj). Consider the continuum limits 
of the sums 

oo 

E i^Ukgy" (51) 

n,fc=0 

and 

oo 

E {A^)n,kgV . (52) 

n,k=0 

The universal behaviors of these two sums give the Laplace transforms of the two functions 
f{n, k) and g{n, k) with respect to z and t. For example 

oo ^ ^ 

E (^)n,fe?'5" dk dn e-^'^'e-^'^'fin, k) = F„(i, z) , (53) 

n,k=0 ^ 

where the subscript "w" indicates the universal part. In order to recover the functions / 
and it suffices to perform an inverse Laplace transform on the universal parts F^it, z) and 
Gu{ti z) of the sums to obtain the functions /(a, /), for which 

F^{t,z) = JdlJ dae-'^-'^fiaJ) (54) 

and similarly g{a,l). Thus 

= . (55) 

f{a,l) 

We shall see this explicitly in the following sections. 

6 Boundary magnetization 

In this section we will apply the discussion above to compute the one- and two-point bound- 
ary magnetizations in the presence of a boundary field. 
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6.1 One-point boundary magnetization 

The boundary magnetization (m) is given by the large n, k hmit of 

(Tr {e'^U - e-^V){e^U + e-^VY-^)n 



{m)n,k 



(56) 



(Tr (e'^f/ + e->'V)>')n 

We may follow the route described in the previous section to compute (m). The first step 



the critical expansion of (j) in the continuum limit, was given in (|37|) . The other quantity we 
need to expand is 

oo 

^, = g0, = ^(Tri?Q'=)g^+^ . (57) 



fc=0 



When h = 0, (pr vanishes by symmetry. When h ^ 0, we can compute (f)r{h) by solving a 
linear combination of the first two loop equations of (plQl). From these it follows that 



(1 



^2h 



c-e^^ + ce^" + ce^'^)(l - 0) + 2pie^^g 



^2h 



4h\ 



g^(l + e"^^ + e 



-,4/1 ^ 



cq + e^^q + e^^q + ce^^q - 2e^^g 
and hence, expanding and multiplying by qcG'"^, we obtain 



(5^ 



(e^'^ - l) (-1 + (3 - V7)e2^ - (4 - 3V7)e^'' 



^e^h _ 1^ + ^1 + 2^7)62'^) 



10 e' 



6h 



2h 



(^e^h _ + (2 + V7)e 

5 • 27/3e2/i (1 + e2ft) 
2^/3 (a + (2 + v^)e2M (l + 4e2^ + e^'^ 



(59) 



1562^^(1 + e^^^)^ 



^5/3 ^ (2T + ^2) $ 

$2 _ (4T)4/3 



+ 0(6^ 



The universal part of ipr is equal to that of 0, up to an /i— dependent constant. There is 
therefore no need to explicitly compute ipr{ci,i), the inverse Laplace transform of ipriZjT), 
in order to determine the boundary magnetization. We need only compute the ratio of the 
universal parts to get the / and A independent result (for h > 0) 



(m) 



_ (e^'- - l) (3 + (2 + v/7)e2^) 

+ i-1 + V7)e^'' + (2 + ^/7)e^^) 



(60) 



Note that (m) is independent of / and A, and is continuous at h = 0. (Note also that the 
expression given in B contains a typographical error in the numerator.) 
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<m> 




Figure 3: Boundary magnetization (m) as a function of boundary field h in fiat space (dotted 
line) and on a random surface (bold line) 

In this particular case, there happens to be a simple argument that gives (m) more 
directly than the computation outlined above. In the large k limit 

(Q')n~gc(/^)-'^7r/(n,A:) . (61) 

Differentiating both sides with respect to h, we obtain 



k{Q'-'R)r^ ~ k [-qcih)-^''-'^qcih)fin, k) + 0{l/k)\ g;^ , (62) 
from which it follows directly that in the large k limit, 

^ ' q,(h) l + (-l + v^)ea + (2 + /7)ei'" ' * ' 

where we have used (PT|) for qdh). This confirms the result obtained in (|5(]|). 

A graph of the boundary magnetization is shown in Fig. 1 (bold curve). As expected, 
with no field the magnetization is zero, and for an infinite field the magnetization is 1. This 
result is compared with the boundary magnetization on a half-plane in fiat space, computed 
by McCoy and Wu (dashed curve). Whereas in fiat space the magnetization scales 
as hlnh for small h, leading to a divergence in the magnetic susceptibility at the critical 
temperature, on a random surface we find that the magnetization is linear at h = 0, with a 
finite susceptibility 

X = 5.Wk=o = ^^|^. (64) 
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6.2 Two-point boundary magnetization 

Having computed the magnetization at a single point on the boundary of the disk in the 
presence of a boundary magnetic field, we would now like to compute the correlation between 
two spins on the boundary of the disk, which are separated by k and / edges in the two 
directions around the boundary. To compute the two point magnetization, 

2^ _ (Tr {e'^U - e~^V){e^U + e-^Vf{e^U - e~^V){e^U + e-^\/)% 
('^ ) - (Tr {e^^U + e-Wf+'+^)n ' ^ ' 

we need expressions for 

q2) = qiq2(^{qi, ^2) (66) 

where 

00 00 

^{qi, q2) = E E(Tr RQ'RQ')q'A , (67) 

fc=0 1=0 

and for 

oo oo 

P(?i, q2) = E E(Tr g^+'+^)gr^g^+^ . (68) 

k=Q 1=0 

Again, using the loop equation techniques of , one can derive the following set of equations 
for a (here, a derivative Dq- denotes a combinatorial derivative as in (|l3), where all variables 
other than are held constant): 

(^{quq2) = &Dg^0(gi,g2) + e(Dg,(T(gi,g2) + Dy,(T(gi,g2)) + tt0(gi)0(g2) 

+70(^1, q2) (0(gi)gi + 0(^2)^2) + bar{qu q2) 

0(gi, q2) = bDl^(f){q2) + 2eDq^(f)r{q2) + &0rr(g2) + aDg^0(gi, g2)gi 
+d (Dg,a(gi, ga) + Dg,(T(gi, gs)) qi + 70(?i)0(?2)gi 

+70(g2)^g2 + /50(gi, g2)gi (0(gi)gi + 0(g2)g2) + agi(Xr(gi, 52) 

0(gi) = 1 + aD5^0(gi)gi + 2dDg^0r(gi)gi + a0„.(gi)gi + (3(t){qifqi^ 
0(52) = 1 + aD^20(g2)g2 + 2dDq^(t)r{q2)q2 + a0rr(g2)g2 + (3<P{q2fq2^ 
(priqi) = &Dg^0(gi) + 2eDg,0,.(gi) + 60„'(gi) + 70(9i)^9i 

0r(g2) = feDq20(g2) + 2eDg20r(g2) + &0rr(g2) + 70(?2)^g2 , 



where 



^^r(gi,g2) = EE(Tri?g'i?g'i?)gtg2 (69) 

fc=0 z=o 
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and 



0(?i,g2)-EE(Trg'i?gV?2- 

fc=0 1=0 



It is easily verified that 0(^1, ^2) is given by 



0(51, 52) 



Qi - q2 



(70) 



(71) 



An expression for a{qi, q-^) in terms of (f){qi), 0(^2) and pi and can be obtained by solving 
these equations. The expression is too long to be included here. On the other hand p{qi, ^2) 
can be directly expressed in terms of (f){qi) and 0(52) in a very simple way as 



. X q2qlBl<P{qi) - qiqm,<P{q2) . ^ q2<P{qi) - gi0(g2) 

p{qi, q2) = = 1 + 



qi - q2 



qi -q2 



(72) 



Armed with expressions for cr(gi,g2) and p{qi,q2), it is then straightforward, if rather 
tedious, to obtain the critical expansions of E and p. They are given by 



E(Zi,Z2,T) =(7e + 



(e2'^-l)'(3 + (2 + v^)e2^ 



$(Zi,T)-$(Z2,T) 
Zi — Z2 



and 



p(Zi,Z2,r) = 



20 • 2V3e4/^ (1 + e^'f 

^z,,T)-^z2,T) 



eV3 + C)(e2/3) 
(73) 



1 



20-2V3a(/?,)2 



^1/3 ^ 0(^2/3) _ 



2/3 N 



(74) 



Zi — Z2 

As in the case of the one-point magnetization, the universal parts of E and p depend on 
Zi, Z2, and T in the same way. Consequently, the ratio of the Laplace transforms of these 
universal parts will simply be the /i-dependcnt ratio of the universal parts themselves. It 
follows that the two-point boundary magnetization. 



-,2h 



1) (3 + {2 + V7)e 



2h 



P (1 + (-1 + /7)e2'^ + (2 + ^/7)e^^) 



2 ' 



(75) 



is precisely the square of the one-point magnetization. We find an absence of polynomial 
corrections to this result (at least to the first subleading order), indicating that correlations 
between boundary spin operators decay exponentially, as one would expect in analogy with 
the flat space theory. 
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7 Bulk magnetization 



The expectation value of the bulk magnetization in the presence of a boundary magnetic 
field, on a disk with boundary length k and area n, is given by 

This can be evaluated by considering cylinder amplitudes with one boundary having a bound- 
ary magnetic field, and the other with a single boundary edge. The second boundary rep- 
resents a marked point on the bulk. Although the second boundary corresponds to only a 
single edge rather than 3 edges as would be appropriate for a triangle corresponding to a 
single spin, this distinction should not be relevant in the continuum limit where the bound- 



ary becomes pointlike. Again, a quantity such as ([76| ) can be computed by the method of 
loop equations 0. 

To compute the magnetization, we require two punctured-disk amplitudes: 

oo 

T{h) = J2 (Tr (e'^U + e-^Vf ■ Tr (f/ - V))q^ (77) 

fc=0 

and 

oo 

\{K) = (Tr {e^U + t^^'Vf • Tr (f/ + V))q^ . (78) 

fc=0 

As in (0)) we define functions related to r but with additional words corresponding to 
sequences of spins on the outer boundary: 

oo 

^«.(.,r)(g) = E(Tr ^(Q, i?)Q' ■ Tr (f/ - l^))g^ . (79) 

fc=0 

The first step in computing the bulk magnetization is now to derive a set of eight inde- 
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pendent equations which close on the quantities (r, r^, r^r, r^gr, ^rrr, Trqqr, ^rqrr, T„rr)'- 
T = 2/5g^0r + aqTrr + aD^r + 2(ir,. + Cgqcp 

Tr = 2jq(f)T + bTrr + bD'^T + 2eDgTr + Cr(f) 

Trr = ar + jq {Tcf)^ + (pTr) + CT^qr + hTrrr + ^D^r^ + eDgT^r + Cr.</)r 

DqTr = 7r + /5g (r^^ + ^r,.) + dT^qr + Cfi^rrr + aD^Tr + dDqTrr + Cq(f)r 

^l^r = 7 i.Pq^ + iq<P) + P^r + /^^ (t^^ + M + 



Trqr = a {pqT + tqCj)) + -fTr + 7g (rDg^^ + ^DgT^) + 
Tj-rr = a {prT + + OT^ + 7g (r0,..,. + (f)Trr) + 



(80) 



where 



tq = {TtQ -Tt {U -V)), tr = {Tt R-Tt {U -V)) . (81) 
As before, Pq = (Tr Q), p^ = (Tr R), and the various constants take the same values as in 



The equations ( ^O]) can now be used to express r(/i) as a polynomial function of (pih). 
However, this equation will also contain a number of unknown correlation functions tq, tqq, 
some of which appear explicitly in (|80|) and some of which arise from 
the derivatives of r, since 

D„r = g-i(r- 1) , 

' ^ ^ (82) 

D^r = g-2(r-l-t^g) , 

and so on. As in the computation of 0, these correlation functions can be reduced to a much 
smaller number of unknowns by expanding (PP]) order-by-order. It turns out that after using 
all the relations in (I5DD (c/. (pi])), two extra relations are required between 

t„ = (Trf/-Trf/), = (Tr f/^ ■ Tr f/) , 

(83) 

U = (Tr ■ Tr U), = (Tr ■ Tr U). 
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The first extra relation comes from tlie calculation in of the critical expansion of 

Wo = ^(Tr (U-V)- Tr {U - V)) = - , (84) 

which is given by 

yj^ = ^^1^ (l + 5'/h'/H' + 0(eV)) . (85) 

The second relation is obtained by differentiating the matrix integral expression for pi with 
respect to g: 

dgPi = . (86) 

Then 

d,P. = (87) 

can be used to obtain the critical expansion of dgPi. 

Armed with these extra relations, we have all the information required to expand r in e. 
We obtain 

^ ^ 2^/^^e(l + 2v/7)^(l + (-1 + v/7)e^^ + (2 + v/7)e^^)Z ,3 

15e2^(l + e2^)(<l>(Z,T) + (4T)2/3) ^ j . yoo) 

As discussed in Sec. H, the quantity of interest is the inverse Laplace transform of the universal 
part of r, given in this case by 

f (L, A) = (l + 2v^)^^?c ^2/3^-5/3^-LVA _ (gg) 

Here, we have introduced rescaled area and boundary length parameters as in (^). 

It is much easier to compute the critical expansion of X{h) since it is directly related to 
4>{h) via 

X{h) = dg(j){h) . (90) 

This gives an expansion 

, 2^/^(-l + 2V7){1 + (-1 + V7)e'^ + (2 + v/7)e^^)A(Z,T) „,,3 _,/3 

81^,6^^(1 + 6^^^) ' ^' ^^^^ 

where 

The leading term in A has an inverse Laplace transform 

A(L, A) = ^-^^^^L'^'A-'^'^''^^^ ■ (93) 
24 



The bulk magnetization in the continuum hmit is then given by 

(M) = ffi^ = . (94) 

A(L, A) 

The numerical values of the coefficients of f and A have exactly cancelled. Although it 
may appear that the magnetization can be greater than one, this formula is valid in the 
continuum limit, for which A ~ ^ 1. We also notice that this form of the magnetization 
is independent of h except for the dependence on the scaling factor a{h) incorporated in L. 
At h = 0, this magnetization is discontinuous and vanishes. One nice feature of this result is 
that it correctly reproduces the scaling behavior expected after the magnetization operator 
has been gravitationally dressed according to the KPZ/DDK description of Liouville theory 



21| , |22| . The gravitationally dressed scaling dimension of the bulk magnetization field is 
A = 1/6. By analogy with the flat space theory we expect that the bulk magnetization 
should scale as (M) ~ c?"^^ where c? is a measure of the distance from the boundary. This 
is precisely the behavior seen in (CT), since A/L has dimensions of length. 



8 Discussion 



Let us summarize the implications of the results we have derived in the previous sections. 
8.1 Renormalization group flow 

In Section 4 we computed the disk amplitude in the presence of a boundary magnetic field as 
a function of the disk area a, the boundary length /, and the boundary field h. We discovered 
that the result could be written in terms of the two variables A = a/5 and L = a{h)l, in 
which case the disk amplitude took on precisely the form of the analogous function when 
the boundary conditions are conformally invariant. Thus, the effect of a boundary field on 
this amplitude amounts to a rescaling of the boundary length by an /i-dependent factor. In 
Section 7, meanwhile, we found that the bulk magnetization as a function of a, I and h could 
also be expressed in terms of A and L, and that its functional form was the same as in the 
presence of an infinite boundary field. 

Taken together, these results imply the existence of an RG flow to the conformally in- 
variant boundary condition with an infinite boundary field; h is a relevant operator which 
goes to ±oo in the infrared (in this context, as the disk area and length grow large). Further 
evidence is provided by the discontinuity in the rescaling function at h = 0: any imposed 



25 



boundary field, no matter how small, leads to magnetization in the bulk of the same form 
as that expected in the presence of fixed boundary conditions. 

A related phenomenon has previously been derived for the fiat-space Ising model on a 
half-plane geometry [0, ||, ^. Again, one can compute the magnetization of a point in the 
bulk in the presence of a boundary field; however, rather than depending on the area of 
the surface and length of the boundary (both of which are infinite for the half-plane), the 
magnetization is a function of the distance from the boundary. (In quantum gravity, where 
we sum over all geometries, it would be conceivable but much more difficult to compute any 
quantity as a function of, say, minimum geodesic distance from the boundary. Computing 
anything "at a fixed point" is even more problematic, and not really well-defined in the 
absence of additional fields.) Chatterjee and Zamolodchikov |^ show that the asymptotic 
form of the bulk magnetization depends on the distance from the boundary as y~^^^ in the 
presence of any nonzero boundary field, just as it does for fixed boundary conditions [p3| . 
Our results demonstrate that this RG fiow is preserved in an appropriate form after coupling 
the theory to quantum gravity. 

8.2 The dual picture 

As discussed in Section 2, the Ising model on a random lattice can also be formulated in the 
dual picture, in terms of matrices X and Y , where X denotes an edge separating two equal 
spins and Y an edge separating two opposite spins. We have seen in the introduction that 
the change of variables 

' (95) 

in the action for the Ising model leads to the dual action. Thus any calculations in the 
original model can be reinterpreted in terms of the dual model. We shall now discuss the 
duals of our results for the boundary and bulk magnetizations in the presence of a boundary 
magnetic field. 

First, let us discuss the boundary conditions corresponding to the weights 

(Tr Q") = (Tr {e^{X + Y) + e-^{X - Y)Y) = (Tr (cosh(/i) X + smh{h) y)") (96) 

where we have dropped factors of for notational simplicity. In the dual variables, h = 
corresponds to fixed boundary conditions, while h = ±oo corresponds to the two types of 
free boundary conditions, {X ± F)". Thus, in this picture h plays the role of a "boundary 
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freeb.c's U+V X+Y, X-Y 

dualities 

fixed b.c's U, V X 

Ising dual Ising 

Figure 4: The duality map that interchanges free and fixed boundary conditions. Note that 
for the Ising model, there are two fixed states and only one free state, whereas the opposite 
is true of the dual model. 

freedom field" rather than a boundary magnetic field. In the limit h ^ the boundary 
condition is fixed (all X's), while in the limits h —>■ ±00, the boundary conditions are a pair 
of free boundary conditions with different signs in the weights assigned to configurations 
with an odd number of F's on the boundary. These fixed and free boundary conditions 
in the dual model are of course precisely the Kramers- Wannier duals of the free and fixed 



boundary conditions of the spin representation. Correspondingly a spin operator in the 
original variables is transformed into a disorder operator in the dual variables. 

Our results in Sec. |^ for the boundary magnetization can thus be reinterpreted as a 
calculation of the expectation value of the boundary disorder operator. We see that as 
the boundary freedom field is increased, so the expectation value of the boundary disorder 
operator 



{Tt {e^{X + Y) -e-''{X -Y)){e''{X + Y) +e-''{X -Y))'') 
^ ^"'^ ~ (Tr {e'^iX + Y) + e-\X - 

(Tr (sinh(/i) X + cosh(/i) Y) (cosh(/i) X + sinh(/i) Y f)^, 
(Tr (cosh(/;.) X + sinh(/i) Y f^^)n 



(97) 



tends to 1. At /i = the boundary condition is fixed and so, by symmetry, the boundary 
disorder must vanish. 

The dualization of the calculation of the bulk magnetization goes along similar lines, and 
we conclude that the bulk disorder in the presence of a boundary freedom field is given by 
the expression 

(D) = L^I^A-^!^ (98) 

except when d = 0, when the bulk disorder also vanishes identically. There is a striking 
consequence of this result. In the Ising model in ?7, V variables, the coupling of a boundary 
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magnetic field was seen to induce renormalization group flow from free to fixed boundary 
conditions. On the other hand, when a boundary freedom field is coupled in the dual 
formulation of the model, the renormalization group flow is from fixed to free boundary 
conditions. This change in the direction of the renormalization group flow is a natural 
consequence of the Kramers- Wannier duality of the system. The duality symmetry implies 
that the ground state degeneracies of the free and fixed states are swapped under the duality 
transformation, and so the reversal of RG flow is consistent with Affleck and Ludwig's g- 
theorem 0]. In general, one expects that the RG flow of the theory should be towards the 
conformal boundary condition with smaller degeneracy, so it is natural that the direction of 
the flow switches under the duality transformation. As we see, this result seems to hold in 
the theory equally well after coupling to quantum gravity. 



8.3 The effects of gravity 

On any fixed lattice, the introduction of an external magnetic field on the boundary leads 
to a direct effect on the Ising spins, as there are no other degrees of freedom with which to 
interact. It is therefore natural to expect that coupling to quantum gravity, which introduces 
the local geometry as an additional degree of freedom, will lead to quantitative changes in 
the response of the spins to the boundary field, and indeed this is what we have observed. 

The calculation of the boundary magnetization in Section 6 can be compared with the 
results that have been obtained in flat space by McCoy and Wu 0. They find that the 
magnetization scales as hlnh for small h, and as a result, the magnetic susceptibility x 
diverges at the critical temperature. On the other hand, we have seen in ( |6^ ) that the 
magnetic susceptibility at the critical temperature is finite when the Ising model is defined 
on a random lattice. It seems likely that the exact numerical value (Q) of the magnetic 
susceptibility x depends on the discretization scheme we have used, and is not universal. 
However, we expect the fact that x is a finite constant to be universal (independent of the 
specific discretization scheme chosen), although in the absence of calculations in alternative 
schemes, this remains a conjecture. 

The effect of gravity is therefore to soften the initial impact of the boundary field. It is 
natural to interpret this softening as being due to the interaction between the spins and the 
geometry; the coupling between spins and the boundary field changes the relative weighting 
of different geometries, which changes in turn the effect of the neighboring spins on any 
one boundary site, leading to a more gradual increase in the boundary magnetization as a 
function of boundary field. 
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(a) (b) 

Figure 5: The relative weight of geometries of different shapes changes with boundary mag- 
netic field h. Roughly speaking, for large h, geometries of type (a) are suppressed relative 
to type (b). 

A related aspect of our results is that the bulk magnetization is seen to decrease when 
a nonzero boundary magnetic field is increased (at least asymptotically, for large areas). At 
first this seems implausible, and indeed at a fixed point in flat space, the bulk magnetization 
cannot behave in this way. However, the sum over geometries provides a possible explanation 
for this unusual effect. The expectation value of a spin in the bulk naturally depends not only 
on the magnitude of the boundary field, but also on the average distance of the point from 
the boundary. Therefore, the decrease in the bulk magnetization can arise if the boundary 
field alters the relative weights of different geometries in such a way as to move a typical 
interior point further away from the boundary (as in Figure 5) . 

We have therefore verified that a number of features of the Ising model in flat space are 
maintained in the presence of quantum gravity, while also demonstrating that the dynamical 
geometry does have a measurable effect. It would be interesting to check more directly that 
the explanations we have given for these phenomena are correct, for example by numerical 



simulation methods such as those described recently in ||2J . 
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